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Abstract. - We propose an explanation for the commonly-seen fading in the density of a stain 
remaining after a droplet has dried on a surface. The density decreases as a power p of the 
distance from the edge. For thin, dilute drops of general shape this power is determined by a flow 
stagnation point in the distant interior of the drop. The power p depends on the local evaporation 
rate J(0) at the stagnation point and the liquid depth h(0) there: p = 1 — 2 (h(0)/h)(J / J(0)) y 
£■ — . where h and J are averages over the drop surface. 
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Introduction. — Recent years have witnessed 
startling new mechanisms for creating predictable, self- 
organized structures from long-known classical principles 
[1]. Among these newly recognized mechanisms is the de- 
position of solute in an evaporating drop of liquid [2]. 
Evaporation entails a specific flow pattern of the liquid 
which in turn concentrates the solute onto the perimeter 
f the drop. Many universal features of the resulting de- 
position have been tested [3, 4] and exploited for applica- 
tions [5-9]. However, several prevalent features of the den- 
sity profile of the deposition still await explanation [10]. 
Here we address one such feature: the decrease or fade- 
out of density with distance from the perimeter, as shown 
in Fig. 1. Our mechanism predicts a power law fadeout 
whose exponent depends on two controllable features of 
the evaporation profile and droplet shape. 

The drop of Fig. 1 contained a volatile solvent that 
partially wet the solid substrate; this solvent leaves the 
drop over time via evaporation. The drop also contained 
nonvolatile molecular or colloidal solutes. These may be 
carried along by any lateral flow of the solvent as it dries, 
but they are not carried away when the solvent evapo- 
rates. Such lateral flow occurs because a) the perimeter 
or contact line is typically pinned by irregularities in the 
surface and previous solute deposition, b) the free surface 
of the liquid takes the equilibrium shape of constant mean 
curvature dictated by its surface tension, and c) the local 
change of volume dictated by the thinning of this equilib- 
rium shape is not matched by the local loss due to evapora- 
tion. To supply the volume needed for evaporation, lateral 
flow is required. Near the perimeter the evaporative loss 




Fig. 1: Optical micrograph of the edge of a 5 mm drop of black 
ink diluted with water, deposited on a glass microscope slide 
and dried in air. Continuous decrease of image density occurs 
as one moves downward from the edge of the drop towards 
the center. The discrete particles and other patterns are not 
addressed here. 



greatly outweighs the supply due to local thinning; thus, 
the lateral flow is strong. It is sufficiently strong to carry 
any point of the drop's interior to the perimeter during 
the drying time [2]. For a thin, circular drop with evap- 
oration controlled by air diffusion, the flow field may be 
readily determined and the consequent accumulation of 
solute with time deduced [3]: in this base case the mass 
deposited at the perimeter in time t varies initially as the 
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While these facts explain the strong accumulation of the 
solute at the perimeter, they do not explain the profile 
of solute density seen in evaporated drops. From a well- 
defined outer edge, this density rises steeply as one pro- 
ceeds inward, reaches a maximum, and then fades gradu- 
ally away. The width w of the deposit increases with the 
initial concentration of solute. To understand this profile 
one must relate the solute's motion to its concentration. 
In 1996 Dupont [11] proposed a simple mechanism to ac- 
count for the effect of concentration. Recognizing that 
solute motion must be inhibited as its concentration in- 
creases, he postulated that the solute simply stops mov- 
ing when its local volume fraction </> exceeds a threshold 
value fa, while solute at lower volume fraction is carried 
passively. Following earlier work by Popov [12], Zheng has 
recently analyzed the surface density profile resulting from 
the Dupont mechanism in two important cases [13]. Zheng 
defines the surface density S as the volume of solvent per 
unit surface area measured after the solvent has evapo- 
rated. For our base case this Y,(x) varies with distance x 
from the perimeter as x~ 7 . 

Below we simplify and generalize Zheng's analysis to a 
broad class of evaporation profiles and droplet shapes. We 
first define and justify the simplified physical system to be 
analyzed. We recall [12,13] that the density profile is con- 
trolled by the advance of the solidification front separating 
the mobile from the immobile solute. We observe that the 
fadeout deposition corresponds to a late stage of evapora- 
tion controlled by a stagnation region of the flow near the 
middle of the drop. We then show that the fadeout expo- 
nent p is entirely determined by a dimensionless "Poisson 
ratio" characterizing the distant stagnation point. We first 
consider circular droplets and then more general shapes. 

A sessile drop wets a circular region of fixed radius X 
contains a solute at initial volume fraction 0, <C <f> c , uni- 
formly dispersed in the fluid. Since the solute is dilute, 
the width w of the deposit will be arbitrarily narrow in 
comparison to the drop radius X. The initial contact 
angle is small, so that the height h(r) at distance r from 
the center is much smaller than X everywhere. We take 
X itself to be small enough that gravitational distortions 
of the shape are negligible. Thus the initial h(r) can be 
written hi{r) = Hi (1 — (r/X) 2 ), where H is the center 
height. At time t = evaporation begins. The volume 
removed per unit time and per unit area from the surface 
at r is denoted J(r) (Fig. 2 inset). This J depends on 
the thermodynamics of evaporation and the environment 
of the drop, but not on time [3]. By increasing the satu- 
ration of the surrounding atmosphere, we may make J as 
small as we wish. We thus choose J to be so small that 
kinetic effects in the liquid such as viscous dissipation are 
negligible: the flow is quasistatic. We choose solute par- 
ticles that do not absorb on the surfaces. Their diffusion 
during the drying time T is made negligible on the scale 
of X but substantial on the scale of H , so that the solute 
remains uniformly dispersed in depth during the evapo- 



ration. Further, the solute particles must be sufficiently 
small or polydispcrse that the solute concentration may 
be treated as a continuum and discrete-particle layering 
effects [14] are avoided. For the moment, we suppose that 
the solute is perfectly compressible, so that the immobi- 
lized solute does not alter the shape of the fluid surface. 

As the evaporation proceeds, the central height H de- 
creases to zero linearly in time (since the rate of volume 
loss is constant). A lateral flow at the depth-averaged 
speed v(r, t) arises in order to restore the imbalance of 
volume loss between J{r) and h(r, t): V ■ (v h) = —h — J. 
Since J(r) and h(r) = —hi(r)/T on the right remain fixed 
in time, the v h on the left side is time-independent as 
well. However, the solute profile far) undergoes qualita- 
tive changes from its initial uniform state as drying pro- 
ceeds. This evolution gives rise to changing deposition 
over time. The two limiting regimes of interest are the 
early regime t/T -C 1 and the late regime (1 — t/T) -C 1. 
In this late time regime the fraction of remaining sol- 
vent is small and solute from the central region will have 
reached the perimeter. On the other hand, the remaining 
solvent may still far exceed the amount of solute. We shall 
suppose that the initial volume fraction fa is so small that 
the average volume fraction remains much smaller than 
4> c during the late as well as the early regimes. We shall 
not consider the final regime, in which (f> has become com- 
parable to <j) c , since it contributes arbitrarily little to the 
deposition profile of interest. Also, by avoiding the final 
time regime we avoid complications like concave meniscus, 
contact line depinning and viscous stresses. 

Deposition front. — If the solute were transported 
purely passively, a nonzero mass of solute would be trans- 
ported to the perimeter over time. Thus the density there 
would immediately become infinite. Imposing the Dupont 
constraint 4> < 4>c prevents this infinite density. Immedi- 
ately a zone of immobilized solute appears at the perimeter 
and begins to expand inward. Outside this immobilization 
front, lateral motion ceases and the surface density 
becomes fixed in time. However, the fixed solute does not 
alter the outward flow profile v h, which depends only on 
the height and evaporation profiles. Under the quasistatic 
conditions we have chosen this v h remains unaltered by 
the solute. 

The motion of the front can be found by a mass-balance 
argument, as sketched in Fig. 2. Denoting its current dis- 
tance from the perimeter by L(t), its advance AL in a 
time At is proportional to the amount of solute carried to 
the front from the interior during At. The solute reach- 
ing the perimeter at time t originated from some radius 
ri(t) in the interior of the drop. In time At, material is 
advected to r-j from a point r\j — v(rj,,t)At just inside r,. 
Thus [3,13] r, = —v(ri, t). The volume AV of solute reach- 
ing the perimeter at time At is thus fa 2nri hi(r{) Ar.i or 
fa 2irrihi(ri)v(ri,t)At. The volume reaching the immobi- 
lization front, just inside the perimeter, is the same in the 
limit L -c X. We determine the motion of L by noting 
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Fig. 2: (color online) Inset: vertical section near the edge of a 
drop. Immobile region (solid color) evaporative current J lat- 
eral velocity v and height h are shown. Main figure: sketch of 
volume fraction <f> vs distance from perimeter x in the Dupont 
model [13]. To the right, the volume fraction is unaffected by 
immobilization. To the left, above the volume fraction (j> c in- 
dicated by a horizontal line, the solute is immobile. Here the 
volume fraction <f>(x)is larger than <j> c only because the height 
at x has decreased after the deposition occurred there. The 
immobilized region is bounded by the immobilization front at 
L(t). The solute volume AV advected to the front in an inter- 
val At is shown by the lower dark-colored bar. This increment 
of solute is sufficient to increase the concentration just ahead 
of the front to <f> c for a width AL indicated by upper light- 
colored bar. Accordingly, the immobilization front advances 
by AL. For the dilute case treated in the text, mobile volume 
fraction is very small compared to <j> c . 



that the volume fraction in the AL increment must be 



AL h(X - L,t) 2ttX = AV 



(1) 



Since h(X-L,t) 
yields 

L L 



2L H/X for L < X , this condition 
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(2) 



The quantities on the right side are known. This equation 
thus determines the motion of the front. Given L(t) we 
may determine the surface density X(a;). This density is 
fixed at the moment t^x) when the front crosses x. At 
that moment £(£(*)) = 4> c h(X - L,t). 

Stagnation flow. — We now consider how X(ce) is 
determined in the late-time limit. In this limit the ini- 
tial radius Ti(t) is close to the center where r = 0. Here 
the fluid's vector velocity v{r) cx r. We may determine 
v(r) by considering the evolution of a small concentric 
cylinder of liquid of radius r and area A. The expansion 
AA in this area is driven by the local evaporation cur- 
rent J(0) and height h(0) (Fig. 3, left inset). In time 
At, the cylinder loses volume Afij = J(0) AAt by evap- 
oration. It loses volume through the decrease in height 
by an amount Afi/j = —A Ah. The loss due to evap- 
oration is not sufficient to reduce the height by the re- 



quired amount. Thus the cylinder's area must expand 
by an amount AA such that h(0) AA = Afl h - Aflj = 
A (-A/i(0))(l + J(0)/A(0)). Evidently h(0) and J(0) are 
related, since J is what causes h to change. This rela- 
tionship can be expressed in terms of the area average 
of J, denoted by J, and the corresponding height aver- 
age h: J = -fi = hi/T. Also h(Q) is given by -h t (0)/T 
so that J(0)/h(0) = -(J(0)/J) (hi/hi(0)). Noting that 
hi/hi(0) = h/h(0), we infer 



AA/A=[l-a] (-Ah/h), 



(3) 



where a = (J(0)/J) (h/h(0)). The relative gain of area 
is proportional to the relative loss of height, as in the 
compression of an elastic solid [15] with Poisson ratio v — 
[1 — a]/2. For typical evaporating drops J(0) is positive, 
so that 2v < 1. In order to have an outward evaporating 
flow with a positive AA, v must be greater than zero. 

This "Poisson ratio" v describing the stagnation flow 
controls the motion of the immobilization front and thence 
the deposit profile S(x). Evidently A = —2vA h/h. The 
—h/h is simply the remaining time T — t = u. From A we 
determine fj via A4/A4 — 2vi(ri)/ri — —2fi/ri. Thus Eq. 
3 becomes — 2f,/rj = 2vju, and cx u" . We may now 
evaluate the front equation (2) noting that /ij(rj) — > Hi, 
H(t) = (u/T) Hi and v{r h t) = dn/du cx u"- 1 . 



L L oc u 2y ~ 



(4) 



so that L cx it" -1 / 2 . Since v < | as noted above, L 
diverges as a fractional power of the remaining time. This 
means that the front extends inward to distances much 
greater than the small deposit width w; our theory ignores 
the final regime where L grows comparable to the radius 
X. 

From the front's motion the fadeout profile follows di- 
rectly. As noted above S(x) = <j> c h(x,tL{x)), so that 
E(ar) - ai^+s)/^-!). This amounts to S(x) - a: 1 " 2 /". 
For the base case described above, h(0) — 2h and J(0) = 
2 J [3], so that a = 1/4 and S ~ x~ 7 , in agreement with 
the detailed calculation of Ref. [13]. For uniform evapora- 
tion with J(0) = J, we find a = 1/2 and X ~ x~ 3 , again in 
agreement with Ref. [13]. Similarly, for a gravity-flattened 
drop with h(0) = h with J(0) = \J, we find a = 1/2 so 
that again X ~ x~ 3 . Finally, a flattened drop with uni- 
form evaporation yields nominally a — 1 and X ~ 
However this case lies marginally beyond the scope of our 
assumptions, since there is no reason here for a global out- 
ward flow of solvent. 

General shapes. — The reasoning above may readily 
be generalized to non-circular drops. As before, the evap- 
orative current J if) is constant in time, and the height 
profile retains its initial shape during evaporation so that 
h(f,t) — (u/T)hi(r). As before, each point at position 
y on the perimeter accumulates solute and has an immo- 
bilization front at L(y,t). We denote the height at the 
front as h(y,L,t). The front advances because of solute 
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Fig. 3: (color online) Sketch of the stagnation flow in an evap- 
orating drop of generic shape shown in the right inset. A seg- 
ment Ay at the boundary is marked in color; its initial pre- 
image at time t earlier is shown by the colored segment at the 
left of the main figure. The wedge-shaped area traced out by 
this segment is shown as a shaded region. Selected streamlines 
from a small circle around the stagnation point are shown as 
rows of dots. The two straight lines at lower right show the 
principal flow directions. Opposite to each straight line is a 
pair of lines, showing the convergence of streamlines in the 
maximum-flow direction and divergence in the minimum-flow 
direction. Left inset shows the advance of a small cylinder of 
fluid (unshaded) to the shaded region at the stagnation point. 
The evaporation current J(0) and height h(0) are indicated. 

advected from the interior. As before each point on the 
perimeter receives solute at time t from a point initially 
at Ti{y,t) in the interior. Again this interior point moves 
inward in time according to fi(y,t) — —v(ri,t) . The ve- 
locity is driven by a scalar pressure field which decreases 
along streamlines. 

The Ti(t) extends inward with time and eventually 
reaches a stagnation point of highest pressure. As with 
the circular drop, this stagnation region is the origin of 
the material that creates the fadeout profile. Here the 
stagnation flow is anisotropic as shown in Fig. 3. 

To determine the fadeout profile, we consider the pre- 
image of a segment Ay at time t, denoted Ay^t). We 
may determine its motion by considering a wedge-shaped 
region of the drop traced out by Aj/i(t) as t — > T (Fig. 
3). We then analyze the increase in the area A A of this 
wedge with time. This area changes for the same rea- 
sons as considered above, so that it too obeys AA/A = 
[1 - a] (-Afe(0)/fc(0)) with a = (J(0)/J)(h/h{0)). The 
solute in the area AA is evidently carried to the bound- 
ary segment Ay in time At. As before, this newly added 
solute creates an immobile region of volume fraction <f> c 
and width AL given by 4> c Ay AL h(y, L, t) = (fnhi(0)AA. 
Since AA is given by the same expression as in the circular 
case, the growth of L with time is also the same: L(t) ~ 
U v ~i (cf. Eq. 4). Finally T,(x) = <p c h(x,y,t) -> x 1 ^/ 01 , 
where a = (J(0)/J){h/h(0)). 

Discussion. — We infer that the fadeout-profilc of an 
evaporating drop is robust and general. Of the many 
known forms of singular pattern formation, the fadeout 



mechanism described here is notable for its remotely con- 
trolled quality. The flow properties at a single point of 
the drop control the power-law profile of deposition in the 
remote region at the perimeter. Under the conditions we 
have assumed, the fadeout follows the same power law over 
the whole perimeter (though the amplitude depends on 
position). Though the flow leading to the fadeout is gen- 
erally difficult to calculate, the predicted power depends 
only on two simple features of the system: the normalized 
evaporation current and normalized height at the stag- 
nation point. These may be deliberately controlled. For 
example, the evaporation profile J(r) can be controlled by 
putting a perforated lid above the drop [8] . Thus one may 
shape materials on small scales where explicit molding or 
machining are not feasible, Another aspect of capillary 
shaping of solutes was recently demonstrated by Vakarel- 
ski et al. [17]. The fadeout mechanism described here im- 
plements a form of shaping that complements Ref. [17], for 
use when smooth and controllable gradients are desired. 

Naturally the accuracy of these asymptotic predictions 
is limited in practice. If one ventures too far from the cir- 
cular shapes and J profiles treated by Zheng, the asymp- 
totic regime likely shrinks. Ultimately the qualitative sit- 
uation of a smooth, outward flow from a single stagnation 
point would be expected to break down. Moreover, in or- 
der to find the controlling ratios J(0)/J and h(0)/h one 
must know the position of the stagnation point. For gen- 
eral geometries, this point is not known a priori. 

The abrupt Dupont rule for immobilization appears un- 
realistic at first sight. In reality, the loss of mobility with 
increasing volume fraction is surely more gradual, occur- 
ring over perhaps an order of magnitude of volume frac- 
tion leading to <f> c . However, in the dilute limit analyzed 
here, this abruptness has little effect. Over the time range 
of interest, the volume fraction throughout the mobile re- 
gion is indefinitely smaller than (f> c ; thus the concentration 
range in question occurs over an arbitrarily small spatial 
distance. Therefore a smoother onset of immobility would 
not have affected the deposition. 

Apart from these geometric effects, one expects mate- 
rial effects to alter the deposition. Brownian motion of the 
solute [5] , adsorption on the surfaces, stratification of the 
solute [8] and depinning of the perimeter [3] must distort 
the deposition to some degree. The final regime of evap- 
oration, where the mobile volume fraction approaches <f> c , 
must also cause distortion. In many experimental condi- 
tions these competing effects can easily dominate the ones 
considered here [5] . Yet the competing effects can readily 
be reduced in order to isolate the essential effects treated 
here. 

Our supposition of a compressible solute requires spe- 
cial discussion. We assumed that the solute is perfectly 
compressible even in the immobilized state. Typical poly- 
meric solutes come close to this ideal, since they interfere 
with each other's motion even when their volume fraction 
is only a percent or less [16]. The polymers can thus resist 
the advective flow v while still compressing vertically to 
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many times their immobilized concentration. However, 
many solutes are not compressible. These would build 
up a deposit whose height could not change with time, 
thus forming a dam around the enclosed fluid. The dam 
perturbs the height profile of the mobile fluid within, con- 
tradicting our assumption [12]. In principle, this altered 
height profile creates an altered flow profile and this com- 
plicates the analysis. It must have a strong effect in the 
final time regime, where the fluid surface would eventually 
become concave. If we avoid the final time regime, then 
the height is only appreciably altered close to the deposit, 
i.e., in a vanishingly small fraction of the drop. We ex- 
pect negligible alteration of the stagnation flow. Thus it 
is plausible that the altered height does not invalidate our 
predictions. 

Conclusion. — The fadeout mechanism described 
here represents a novel way to shape material in space 
using the generic properties of stagnation flow in an evapo- 
rating drop far from the site of deposition. The mechanism 
is broadly applicable to a range of evaporation geometries 
to produce a controllable range of power law deposition 
profiles. Further effects eg. from solute dispersion and 
packing will no doubt enrich the possibilities. Such ef- 
fects can be calculated explicitly and offer the prospect of 
further control. 
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